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flj Abstract: Motivated by recent work on the utility of Mellin space for representing conformal 

i-S^ correlators in AdS /C¥T, we study its suitability for representing dual conformal integrals of 

the type which appear in perturbative scattering amplitudes in super- Yang-Mills theory. We 

►^ discuss Feynman-like rules for writing Mellin amplitudes for a large class of integrals in any 

CNj dimension, and find explicit representations for several familiar toy integrals. However we 



show that the power of Mellin space is that it provides simple representations even for fully 

\^ massive integrals, which except for the single case of the 4-mass box have not yet been 

C^ computed by any available technology. Mellin space is also useful for exhibiting differential 

^:^ relations between various multi-loop integrals, and we show that certain higher-loop integrals 

T— I may be written as integral operators acting on the fully massive scalar n-gon in n dimensions, 

^ whose Mellin amplitude is exactly 1. Our chief example is a very simple formula expressing 

k>( the 6-mass double box as a single integral of the 6-mass scalar hexagon in 6 dimensions. 
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1 Introduction 

The interplay between conceptual and technical advances has always been an important 
catalyst for progress in theoretical physics. A spectacular example of this has been the 
development of our understanding of AA = 4 supersymmetric Yang-Mills (SYM) theory [1] — 
the exemplar of four-dimensional quantum field theories — especially since the discovery of 
the AdS/CFT correspondence [2] almost 15 years ago. In particular, two related aspects 
of planar SYM theory which have generated considerable attention within the past several 
years include the application of powerful integrability techniques [3] for determining operator 
dimensions and the discovery of remarkable mathematical structure in perturbative scattering 
amplitudes [4]. 

More recently some attention has focused on the important problem of understanding 
better the structure of correlation functions in SYM theory, which at strong coupling may be 
computed via AdS/CFT using Witten diagrams [5]. Sometimes, both conceptual and tech- 
nical progress can be aided by 'using the right language' (or set of variables), as dramatically 
evidenced for example by the use of spinor helicity and momentum twistor variables for scat- 
tering amplitudes. Motivated in part by appreciation of this lesson, it has been suggested [6] 
that the 'right' place to study AdS/CFT correlation functions is not in position space but 
rather in Mellin space, the benefits of which for correlation functions in general CFTs were 
pointed out in the pioneering work of Mack [7]. 

For the purpose of studying AdS/CFT correlation functions Mellin space is both healthy 
and great tasting — it allows for a dramatic simplification of otherwise intractable computa- 
tions (see for example [8-10]), while at the same time providing a definition for a dual bulk 
S-matrix via the vanishing curvature limit of the AdS/CFT correspondence [6, 8, 11]. Mellin 
amplitudes also make the physics of correlators transparent, by showing in a simple way their 
conformal block decomposition, and their associated OPE coefficients. This is analogous to 
the way in which momentum space clarifies the physics of weakly coupled field theories. It 
seems therefore that a strong case can be made that Mellin amplitudes are the right object 
to consider in any conformally invariant setting. 

Motivated by the promise of this approach, our goal in this paper is to carry out a first 
study of the suitability of Mellin space as a language for the the weak-coupling expansion 
of scattering amplitudes SYM theory. Let us emphasize right away that we are interested 
here in the boundary flat space S-matrix, rather than the bulk S-matrix for supergravity or 
string theory in the flat space limit of AdS. The latter has been studied since the earliest 
days of AdS/CFT (see for example [12-14]) and may be computed as a certain limit of SYM 
theory correlation functions (see [8, 15-17] for recent work). The former may be computed 
at strong coupling via AdS/CFT by introducing a probe brane to provide the gluon degrees 
of freedom [18], but in this paper our interest lies in the perturbative expansion at weak 
coupling. 

Although we consider examples of integrals which may appear in various field theories, 
SYM theory is apparently unique amongst four-dimensional field theories in that in that all 



of its scattering amplitudes are amenable to Mellin representations of the type we discuss 
since they all possess dual conformal invariance [19-21]. This symmetry of SYM theory was 
first noticed in some examples [22] and proven to be a property of the integrand for general 
amplitudes to all loop order in [23]. 

One of our motivations for seeking a new language for loop integrals in flat space is that 
despite remarkable recent advances, actually carrying out multi-loop integrals remains a very 
challenging task for which there is no practical general algorithm (see [24] and references 
therein for some of the most modern magic). This stands in stark contrast to the situation 
for the planar integrand of SYM theory, which is amenable to powerful generalized unitarity 
techniques [25-27] and which can in principle be computed for any desired process via the re- 
cursion described in [28] (see also [29, 30]). Given the current relative simplicity of computing 
integrands but the relative difficulty of computing integrals it is natural to wonder whether 
there exists some kind of 'stepping stone' in between these two quantities. Initially we should 
not necessarily require this stepping stone to exist for arbitrary amplitudes in any random 
field theory, only for the very special amplitudes of planar SYM theory. However we should 
require it to be completely canonical — both the integrand and the integrated amplitudes of 
SYM theory are mathematically well-defined objects which look identical to us today as they 
will to an alien civilization a billion years in the future, and the same should be true of any 
good stepping stone. 

One important example of something halfway between an integrand and its integral, for 
those integrals which can be expressed in terms of a certain class of generalized polylogarithm 
functions, is the 'symbol' described for example in [31-33] and first used for SYM theory 
amplitudes in [34] . Symbol technology has proven very useful in several applications (see for 
example [35-48] ) , but it appears that sufficiently complicated amplitudes even in SYM theory 
involve elliptic functions (of a type familiar in the QCD literature, see for example [49]) which 
are outside the class treatable by current symbol technology. 

In this paper we propose that Mellin space representations of the type recently employed 
for AdS correlation functions might provide useful also for flat space dual conformally in- 
variant amplitudes. After quickly introducing the class of integrals under consideration in 
section 2 and the definition of the Mellin amplitude in section 3, we explore aspects of this 
proposal via several examples in the subsequent sections. In particular, we propose in sec- 
tion 5 that there are simple Feynman-like rules for directly obtaining the Mellin amplitudes 
corresponding to dual conformal integrals with trivial numerator factors. These rules are 
considerably simpler than the ones proposed in AdS/CFT, and we show how the former 
seemingly derive from the latter. In practice they amount to thinking of the position space 
diagram as a diagram in Mellin momentum space; in this way internal lines map to propagator 
factors and contact interactions have Mellin amplitude equal to one. This remarkable fact 
has important consequences for deriving differential equations for dual conformal integrals, 
and also allows one to immediately solve such equations, thereby reducing the computation 
of certain higher-loop diagrams to simple integrals of one-loop n-gon diagrams. 

In sections 6 and 7 we investigate the Mellin representation of conformal integrals with 



numerator factors, and in particular the pure chiral integrals studied extensively in [23]. We 
find several interesting features: firstly, the magic numerator factors appearing in the chiral 
pentagon and hexagon integrals make some of the Mellin integrals "collapse" onto boundary 
poles. We explain how this works in detail in a simple example, and then use this property 
to derive a representation of the chiral pentagon in terms of derivatives of hypergeometric 
functions. Next we show how for the chiral hexagon, double numerator factors translate into 
second-order differential operators in the Mellin representation. Such differential equations 
have been studied in [37, 50], and our results make it clear that they easily generalizes to 
more complicated examples. 

We finish this paper with a short discussion, followed by various appendices containing 
additional technical results and details. 

2 Setting up 

2.1 Dual conformal integrals in SYM theory 

Here we provide a quick introduction to the important features of the integrals which appear 
in SYM theory loop amplitudes. The most important property, which has been proven to 
hold to all orders in perturbation theory for all amplitudes, is that the integrand is invariant^ 
under dual conformal transformations. We remind the reader that dual conformal trans- 
formations [21] are nothing but ordinary conformal transformations on the dual variables 
xi, . . . ,Xn related to the momenta of the n scattering particles by 

Pi ^ Xi Xi-)-l, Xij = [Xi ^j) \^'^) 

where all subscripts are taken mod n. 

We take this opportunity to immediately break with the standard conventions of the 

amplitude community in order to streamline the notation for this paper. Instead of using dual 

variables Xi which implicitly are null separated from their neighbors {i.e. (xj — Xj+i)^ = 0), 

we will only use as many dual variables as external faces in any diagram under consideration, 

and furthermore let them initially take generic values. For example consider the four-mass^ 

integral shown in Fig. 1. In the amplitude convention this diagram would be labeled with 

eight x's, but since the value of the integral only depends on the four x's shown in (a), we can 

economize the notation by relabeling the diagram as shown in (b). The value of this integral 

is then 

r fax 3;;^3X24 ,^ ^. 



J — L j m"^ {x — x\f-(x — X2)^{x — x-i)^{x — x^ 



^Strictly speaking it is only covariant, but it is rendered invariant after dividing by the tree-level MHV 
superamplitude, which we follow standard convention in doing. 

^We remind the reader unfamiliar with this terminology that the label 'four-mass' is used because the sum 
of the external momenta entering each of the four corners of the box is non-null. There is no ~\-rr? in any 
propagator and there is no breaking the conformal (or dual conformal) symmetry of SYM theory. 
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Figure 1. The one- loop 8-point four-mass integral, labeled according to the usual amplitude conven- 
tion in (a) and according to our streamlined notation in (b). In (a) it is implicit in the notation that 
each Xi should be null-separated from its neighbors Xi-i and x^+i. In contrast the Xi in (b) and (c) 
are arbitrary, and (a) is recovered by a simple relabeling. This integral corresponds in the dual (Mellin 
momentum) space to a tree-level contact interaction (c). ^3 X4 
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Figure 2. The two-loop four-mass double box integral (a) is a particular limit of the 'fully massive' 
double box (b), computed in Mellin space as an exchange diagram contribution to a tree-level 6-point 
correlation function (in blue). The integral (a) is recovered from (b) by taking the limit X4 — >■ X3, 



xe 



Xl and then relabeling 2:5 -^ X4. We define the integral (b) to include the overall factor : 



■14-^25-^36 



in order to provide dual conformal invariance. This reduces to the factor x\^X2i for integral (a). 



Note, importantly, the inclusion of the overall factor 2:^3X24 into the definition of this integral. 
This factor, which is required for dual conformal invariance, is the first reminder that the 
integrals under consideration here are not exactly those of scalar (j)'^ theory, though in many 
cases they are very closely related. 

Of course (2.2) is manifestly the same (again, up to the factor x\'^x'^^ integral which 
computes the tree-level position space CFT correlation function {0{xi)0{x2)0{x'i)0{xi)) 
of four operators with dimension A = 1 interacting via a four-point contact interaction, as 
shown in Fig. 1(c). Henceforth we will always draw the dual diagram (c) in blue directly on 
top of the corresponding integral (b) in order to save space. 

In our approach it is most natural to always begin with the fully massive version of any 
integral under consideration, and then to recover other versions of the integral by taking 
appropriate limits. For example, the two-loop four-mass diagram shown in Fig. 2(a) is a 
perfectly nice finite and dual conformal invariant integral, but we represent it as a limit of 





Figure 3. 'Magic' numerator factors are denoted graphically by a red line crossing an internal face. 



the fully massive integral as shown in Fig. 2(b) and given by 
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(2.3) 



The fact that fully massive integrals are often the simplest to work with in Mellin space is 
one of its most attractive features, since it is opposite to the experience of amplitudeologists 
to whom more massive integrals are necessarily more complicated. For example, while it 
remains an open challenge to evaluate the integral in Fig. 2(b) at just two loops (it is believed 
to involve elliptic functions), even the L-loop generalization of the integral in (a) was fully 
evaluated long ago [51] in terms of standard polylogarithm functions. 
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We will see below in section 5.1 that it is trivial to write down the L-loop generalization of 
the fully massive integral in Fig. 2(b) in Mellin space, from which (2.2) would follow as a 
special case. Let us however temper our enthusiasm (slightly) by pointing out that taking 
such limits of interest is often but not always a trivial task in Mellin space, as we discuss 
below in section 4.3. 

The correspondence between the integrals appearing in SYM theory and in (/>^ theory 
only holds for the simplest diagrams. General integrands in SYM theory have non-trivial 
numerator factors. A particularly nice collection of such integrals are those involving chiral 
numerator factors of the type discussed extensively in [28, 52]. Denoted graphically by a 
red line crossing some internal face (see Fig. 2.1), the corresponding numerator factor is 
proportional to'^ the quantity (x— y)^, where y is a solution to the leading singularity equations 



{y - Xif = {y- Xi+if = {y- Xjf = (y- Xj+if 



0. 



(2.6) 



By convention chiral integrals are normalized so their nonzero leading singularity is 1. 
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Figure 4. The chiral pentagon (a) and hexagon (b) integrals under consideration in this paper. 



For given x's these equations have two different solutions for y; the corresponding two nu- 
merator factors were denoted by squiggly and dashed red lines in [28]. The utility of these 
perhaps strange-looking numerator factors is precisely that by killing one of the leading sin- 
gularities associated with the loop integration variable x they allow one to express various 
integrands in SYM theory very compactly (especially for example MHV amplitudes, as those 
are maximally chiral). 

In this paper we will study in detail one-loop chiral integrals with one and two numerator 
factors, two particularly simple examples of which are shown in Fig. 4. Here we have drawn 
the pentagon with two massless corners and the hexagon with all massless corners. We will 
actually begin both cases by considering an arbitrary (fully massive) one-loop n-gon integral 
with one or two numerators and then take appropriate limits to reach these two special cases. 
Our interest in them in particular stems from the fact that these are the configurations in 
which the integrals enter one- loop MHV and NMHV amplitudes in SYM theory [23] . In these 
limits the numerator factors simplify and the integrands can be written rather simply in terms 
of momentum twistors [53]. For the pentagon we need at least 8 legs to provide the 3 massive 
corners. For later use let us choose to label the legs by their momentum twistor variables Zi, 
i = 1, . . . , 8 in this case as 
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while the massless hexagon is simply 
3 4 
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These two integrals were evaluated explicitly in [50] with the results 

= Li2(l - uia) + Li2(l - U35) + Li2(l - uu) 

- Li2(l - U13U35) - Li2(l - U13U14) + log(n35) log(ni4), 




(2.9) 




vr 



Li2(l - Uu) + Li2(l - U25) + Li2(l - use) + log(tt25) log(u36) - y, (2.10) 



where 



^i,j+i^i+i,j 
^i,i^j+i,i+i 



(2.11) 



Since the formalism we will employ works in arbitrary dimension we will also encounter 
higher dimensional integrals, including the six-dimensional scalar hexagon 
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This integral has been evaluated in special cases including zero [36, 37], one [38] and three 
masses, but not yet for completely general Xj (though its symbol is known [40]). The im- 
portance of this integral for SYM theory scattering amplitudes, and its relation to the inte- 
gral (2.8) has been explored and emphasized in [37]. 

2.2 Ambient space formalism 

The calculation of conformally invariant integrals is conveniently performed in the embedding 
or ambient space formalism, which goes back to Dirac [54]. For a more recent reference 
with several details on the four-dimensional case see [55], and for some interesting recent 
applications to CFTs see [56, 57]. The basic idea is that we can parameterize d-dimensional 
Minkowski space by projective light-cone coordinates in d -|- 2 dimensions. The SO{d, 2) 
invariance group in d -|- 2 dimensions is precisely the same as the conformal group in d 
dimensions. In this way conformal transformations on the coordinates of the d-dimensional 
space become simple rotations of the d + 2-dimensional coordinate vectors. 



More concretely, consider null vectors in d + 2 dimensions which we will invariably denote 
by capital letters P, Q, Y and so on. That these are projective null vectors means we have 

P^'^Pm = -P+P' + P^'P^ = 0, P ~ AP (2.13) 

where M, /x are d + 2 and d-dimensional indices respectively. 

To obtain coordinates in d-dimensional Minkowski space we need to define a reference 
vector, call it /, thereby explicitly breaking conformal invariance. Then the vector P^ = 
P^^ /{—P • I) parameterizes d-dimensional flat space. In practice we make the convenient 
choice 

pM 

-^^ = ^/2(l,x^x'^), (2.14) 

where we have used light-cone coordinates for the first two entries of the vector. The above 
choice is equivalent to setting —P-I= ^P+. By contracting two independent P vectors we 
obtain 

In practice we will drop factors Pi ■ I throughout this paper. These can always be recovered 
by demanding that amplitudes should be invariant under Pi — ?■ aPi, but in any case they 
always cancel out of any conformally invariant expression. 

The embedding formalism is of course well-defined in any dimension. However it is 
important to notice that using null vectors for describing positions allows us to use vari- 
ous spinor-helicity formalisms for particular cases. In d = 4 for instance one can use the 
equivalence 50(4,2) ~ 5'C/(4) to rewrite six-dimensional vectors P^^ as bi-twistors pl^-^l, 
where A, B are 5^7(4) valued. Since these vectors are null, the corresponding bitwistors are 
simple — that is, they can be written as a product of two 5C/(4) valued spinors^ 

pM ^ pAB ^ ,<^^z^Z^ = z\^Zp , (2.16) 

which establishes a map from a d-dimensional vector x^ onto a line in twistor space. However 
for most of our calculations it will suffice for us to use the embedding formalism, which has 
the advantage of being applicable in any dimension. 

3 Definition and properties of MelHn ampHtudes 

In a conformal field theory the two-point function is uniquely determined (up to an irrelevant 
normalization constant) to be 

{4>a{PMa[P2)) = ttAa. (3-1) 



These are of course nothing but Hodges' momentum twistors [53]. 



Two-point functions of fields with differing A are identically zero. The parameter A is known 
as the conformal dimension of the field (p. The three-point function is also uniquely determined 
by conformal symmetry up to a coupling constant, 

(0Ai(Pi)0A3(^2)<Aa3(^3)) =Ca,A,A3 n V" (3.2) 

with (5i2 = ^~^ 2~ — ^' ^"^^ cyclic permutations thereof. Higher-point correlation functions are 
determined up to an arbitrary function of cross-ratios, which are homogeneous combinations 
of internal products of Pi vectors. For instance, a conformally invariant four-point function 
can be written in the form 

A1-A2 A3-A4 

(</'A,(Pl)<^A.(P2)<^A3(^3)'^A4(A)) = (-PW^l^) ; (-PWAs)^ ^ ^(^^^^ (3^3) 

(Pl2)^^(P34)^^ 

where u and v are the cross-ratios 

A2-P34 -P14-P23 ,„ -x 

^=Tr7^' ^=7r7^- (3-4) 

Typically the function F(tt, v) is a very complicated object. This is not unexpected, since 
we know that correlation functions in position space do not usually have a simple structure, 
even for weakly coupled field theories. However, in that case we know what we should do: 
instead of working in position space we Fourier transform to momentum space. There the 
analytic properties of the amplitude are very simple, and their physical meaning is clear. At 
tree-level one sees simple poles corresponding to single-particle states, and branch cuts at 
loop level corresponding to multi-particle exchange. For a generic conformal field theory, the 
Fourier transform is not so useful. This is simply because typical CFTs are strongly coupled, 
and as such we do not expect that going into a basis of approximately free momentum 
eigenstates should help. Perhaps the easiest way to see this is to simply consider the Fourier 
transform of the two-point function (3.1). This behaves as ~ p-2(d-A) ^^^^ g^ ^^ g^g that 
even the two-point function already shows branch cuts for generic A. 

The situation therefore might seem hopeless, but the large symmetry of conformal field 
theory comes to our rescue. Following Mack [7] we introduce the Mellin amplitude M{5ij) of 
a conformal correlation function of scalar fields via the definition^ 

(</>Ai(Pi)...0A„(P„)) = jdS,, M{6,j) llTid,,)Pr/^\ (3.5) 

Let us look at this formula in detail. The main object in the above is M(6ij), the Mellin 
amplitude. We conventionally defined it such that an overall product of F functions always 

""We remark that in momentum twistor language, this is a Mellin transform with respect to 4-brackets 
of the form {ii+lj j+l) ex x^j. However, from a twistor theorist's point of view, it might be interesting 
to contemplate instead a Mellin transform with respect to the variables Zi ■ Wj = (ij — 1 Ji+1)- We thank 
D. Skinner for this comment. 



10 



appears, and this will be convenient later. The Mellin amplitude is a function of the complex 
parameters 5ij which are being integrated over a suitable^ contour in the complex plane. 
These parameters are symmetric in their indices, and satisfy the constraints 



-Ai, 5^5y=0. (3.6) 



For more precise details on the measure we refer the reader to appendix A. Overall there are 
^2 independent parameters, and this is precisely the same as the number of independent 
cross-ratios. This is not an accident, as the constraints (3.6) are precisely those of conformal 
invariance. This can be seen for instance by performing an inversion on the x^ vectors, or 
more simply, by demanding that under Pi — )• aiPi the overall amplitude scales like a- *. Upon 
solving the constraints, the Mellin representation becomes simply a product of familiar one- 
dimensional Mellin transforms, one for each cross-ratio. In particular, contour prescriptions 
are exactly the same as those for the one-dimensional transform, and the inverse Mellin 
transform is simply the product of the one-dimensional inverses^. 

A nice way of thinking about these 5ij parameters is as internal products of momenta [7] . 
Indeed, if we parameterize the 5ij as 

Oij = Ki ■ Kj [6. I ) 

then the constraints (3.6) are automatically satisfied if kf = —A, and ^j ki = 0. In this way, 
one can think of a Mellin amplitude as depending on Mandelstam-like variables built out of 
these momenta, each momentum being associated with a given field. The resemblance of the 
6ij to momenta is not an accident and we will comment on it shortly. 

What have we gained by using the representation (3.5)? As it turns out we have gained 
considerably. Firstly, in contrast to the complicated functional dependence of general correla- 
tors in position space, Mellin amplitudes are simple meromorphic functions of their arguments 
Sij . Secondly, the poles are directly related to the conformal block decomposition of the cor- 
relator. Consider for instance the four-point function. Applying the OPE decomposition in 
the (12) channel. Mack [7] has proven that the full Mellin amplitude {i.e. including gamma 
functions) has poles in the Mandelstam variable s = —{ki + A:2)^, 



p,n 



s — (Ap — Ip) — 2n 



The position of the poles corresponds to the twist of the primary operators (labeled by p) 
appearing in the conformal block decomposition of the correlator. The extra n summation in 
the above correspond to contributions from descendants of these fields. Further, the residues 



''In all examples we consider, the integration contours may be taken parallel to the imaginary axis, with 
real parts chosen so that ^{Sij) > for all i,j. 

'^For sufficiently high number of legs there are ambiguities in the Mellin amplitude due to the vanishing of 
Gram determinants. In practice there is no cause for alarm, but see the discussion in [16]. 
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of these poles give products of the three-point coupUngs of the theory Ci2p,nj up to a known 
polynomial in the remaining independent Mandelstam invariant t = {ki + k^)"^ , a polynomial 
of order Ip, the spin of the primary field p. In this way, the Mellin amplitude makes the 
physical content of a given correlation function manifest: by examining its poles and their 
respective residues we can immediately determine which primary operators are involved and 
what their three-point couplings are. 

The Mellin amplitudes defined above are valid for any scalar conformal correlation func- 
tion. In practice, most of the work on this topic has so far focused on correlators computed 
with the help of the AdS /CFT correspondence. There, passage into Mellin space has allowed 
for a complete solution of the computation of tree-level correlation functions in arbitrary 
scalar field theories [6, 8, 9]. One finds that there are Feynman rules for directly evaluating 
the Mellin amplitude given a Witten diagram. These Feynman rules look remarkably similar 
to momentum space Feynman rules, once we set 6ij = ki ■ kj. In this note we will focus on the 
computation of conformal integrals in flat space, but we will find the same structure at work. 
In particular, there seems to be a set of Feynman rules for the Mellin amplitudes of these 
conformal integrals, and the rules are considerably simpler than those found for AdS /CFT. 

The resemblance of Mellin amplitudes to momentum space scattering amplitudes has been 
understood in the context of AdS /CFT. By considering the large energy limit of scattering in 
AdS, it is possible to show that the Mellin amplitude becomes precisely the d + 1-dimensional 
flat-space scattering amplitude [6, 8]. This means that we expect that there should be Feyn- 
man rules for computing the Mellin amplitude directly, and indeed this is the case. The rules 
were written in a compact form in [9] and shown to obey a BCFW-type recursion relation 
in [6] . These rules originally apply to the computation of the Mellin amplitude corresponding 
to tree-level Witten diagrams, but it is perhaps not too far-fetched to expect that similar 
rules should hold for conformal flat space integrals, and indeed we will find this is so. 

4 Conformal integrals in Mellin space 

4.1 The polygon 

We begin our labors with the calculation of the polygon integral in momentum space — a star 
in the dual position space. This integral is given by 

The integration should be understood as Jd Q = f_°°dx (with Q = v2(l,a;^,x^) as 
explained in section 2.2). We have also defined the convenient shorthand 

h^^. (4.2) 
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The standard way of performing such integrals is to introduce Schwinger parameters, one per 
denominator factor, 



l''^ l^if ' / dQ exp (q . Y, Up)j . 



In = 7T-'' I ^tf » / dQ exp ( Q • "^ UPi ] . (4.3) 

The Q integral is gaussian, and we get 

(4.4) 



Now the point is that if the original integral satisfies the conformality condition Y17=i ^« ~ ^' 
we can drop factors of ^ tj. A nice way to see this is to introduce a partition of unity 



r+oo 

1= / dv5{v-y2u). (4.5) 

Jo . 

After a rescaling of the ti we get 



If we now did the v integral we would recover the usual Feynman parameter form of the loop 
integral, with the ti playing the role of Feynman parameters. Instead, we will do another 
rescaling, tj — )• tij ^ and perform the v integral to get 

^- = ^^1 n 17*^ ( E ^A «-p I E ^^oP^, 1 • (4-7) 

We see that if the conformality condition is satisfied there is a simplification, and our integral 
becomes 



/■+00 dt i \ 

i * \i<j J 



(4.8) 
To proceed we use Symanzik's trick [58], which tells us that 

I n 17*f' «^P(- E i^t^P^^) = 2 / d'^'^- n nS^^)P^'' (4.9) 



Kj Kj 



with 6ii = — Aj, ^^ (Jjj = 0, as explained in appendix A. These are the same 6ij that we have 
introduced in section 3. Recall that these constraints are solved by defining Sij = ki ■ kj, with 
kf = — Aj. Going back to our integral, we finally conclude that 

In= <f d6i,llT{6.,,)Pr/^\ (4.10) 
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Comparing this expression with the definition of the Mellin aniphtude (3.5), we come to the 
happy conclusion that the Melhn transform in the case is simply one, M{5ij) = 1. 

The reader might feel somewhat cheated by this result, but despite the simplicity of the 
Mellin amplitude the integral (4.10) is definitely non-trivial; after all, there are still the factors 
of T{5ij) left, which by convention we have left out of the definition of M{5ij). However, this 
convention is useful since all conformal integrals we shall compute will always include the 
very same product of gamma functions. Their presence is related to the existence of "double 
trace" operators in the conformal block decomposition of the correlator, which always appear 
in such weak coupling computations^. This means knowing that the Mellin amplitude is one 
for all such conformal integrals is useful: it tells us how to straightforwardly write down 
their Mellin-Barnes representation without any further thought. In particular, the 4(i box 
integral (2.2) is simply 

P13P24 j d<5y J] T{8ij) Pr/'^ , (4.11) 



and similarly the fully massive n-gon integral in n dimensions is a similarly simple 2^{n — 3)- 
dimensional Mellin integral. In order to connect with previously known results from the 
literature we consider in detail the cases n = 4 and n = 6 in appendices C.l and C.2 
respectively. 

There is another more interesting reason why we have chosen the convention (3.5) for 
defining the Mellin amplitude. This is because the fact that the Mellin amplitude as we've 
defined it is simply a constant (in this case) is analogous to the fact that in momentum space, 
non-derivative contact interactions also correspond to trivial amplitudes — they are constants 
as well (up to momentum-conserving delta functions). The natural question that arises is 
whether this analogy will continue to hold when we have a conformal integral corresponding 
to an exchange diagram. An example of such an integral is the scalar double box integral in 
4 dimensions, which is dual to a position space 3—7-3 exchange diagram in (j)^ theory, shown 
in Fig. 2(b). It is precisely to this kind of integral that we shall consider next, and indeed we 
shall find that this analogy continues to hold perfectly. 

4.2 Two polygons with one common edge 

Let us consider then a two-loop computation. Once again we will keep our calculations 
generic, and consider not a double box but rather an n-gon and an ?7i-gon, as before, glued 
together along one edge. The integral in the dual position space is now 



^"'' f.r^.r, tV r(A,) ^ r(A 

{-P^ ■ Qi)^» ii i-Pj ■ Qi)^^ i-Qi ■ Q2 



^Q^^Q^ n (Z^^ n (Z^^f^TZT^^T^- (^-12) 



where L, R refer to the left and right polygons. The only restriction is that both integrals 
are conformal, so that X^iei ^i + ^ = J2ieR \ + ^ = d. As before we introduce Schwinger 



For a discussion of these issues in AdS see for instance [59]. 
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parameters and compute the Q integrals. Conformality simplifies matters and we end up 
with 

r(d) Jo ^^ u Jo s 

with Pl,r = "^ieL R^i^i- Exactly as in the previous case, we use Symanzik's formula (4.9) 
to directly derive the Mellin representation from the above integral. But this time, because 
of the "internal" Schwinger parameter s, which came from exponentiating the {—Qi ■ Q2)~ 
factor, we get a non-trivial Mellin amplitude 

1 r+co J 

r(o) 7o s 

= 1 r"^/-L(i + ,2)-i(E,e.A.-.,) (4_i4) 

r(o) Jq s 

with sl = — (X^jgi ki)"^. Performing the integral we find the simple result 

rf^)r(^) 

M{sl) = ^ ;/, , . (4.15) 



To better exhibit the pole structure of the Mellin amplitude, we can write it instead as 

n=0 

This result is in precise accord with the general expectations of section 3. As predicted 
by Mack, we find poles corresponding to the primary field of dimension 5 being exchanged, 
plus an infinite series of poles labeled by n corresponding to its descendants. The Mellin 
amplitude therefore looks like precisely a momentum space scattering amplitude, in this case 
an amplitude involving exchange of an infinite set of massive fields, with particular propagator 
normalization factors. 

We should not forget however, that at the end of the day our interest is in analyzing the 
conformal integrals appearing in SYM theory scattering amplitudes. In particular, the fully 
massive double box of Fig. 2 corresponds to setting 5 = 1 and d = A. When we do this, the 
infinite sum beautifully collapses onto the single term 

M{sl) = — ^ (4.17) 

1 -SL 

Why is this so? Notice that the first argument of the Pochhammer symbol {1 + 5 — h)n 
vanishes ii 6 = h — \ = d/2 — 1: in other words if the conformal dimension 5 is that of a free 
field. But for free fields we have 'dcj) = 0, which effectively kills off the contributions from 
descendant fields to the expression (4.15). 
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To summarize, we have found that the fully massive two-loop double box integral, which 
appears to be too difficult to evaluate in position space with current methods, can be computed 
by substituting 

SL = -ih + k2 + hf = 3 - 2{5i2 + Si3 + S23) (4.18) 

into (4.17) and (3.5), yielding the extremely simple representation 



IPuP.n,fds, J^ Ylris,,P;^"' (4.19) 

as a Mellin integral of dimensionality 26(6 — 3) = 9. In section 5.1 we argue that this result 
generalizes to higher loops in a straightforward way to provide a (2L — 1)(L + l)-dimensional 
Mellin representation of the fully massive L-loop ladder integral. 

4.3 Massless limits of Mellin amplitudes 

We have emphasized that for Mellin amplitudes it seems most natural to always begin with 
the most massive version of any integral of interest. This means that all Xi are taken to be 
arbitrary, which from the point of view of conformal correlation functions is the most natural 
thing to do. However many of the integrals appearing in scattering amplitudes require special 
constraints on the kinematics — often one or more pairs of cyclically adjacent x's are lightlike 
separated from each other {xjj^_^_^ = 0)) or even exactly equal to each other (xj = Xj+i). The 
question naturally arises then as to how one can specialize to such cases since the definition 
of the Mellin amplitude (3.5) seems to break down if any x?- = 0. 

In many cases (including all of the examples we consider explicitly in the following sec- 
tions), taking the limit as some xf- — ?■ can be done quite easily at the level of the Mellin 
amplitude. This happens whenever the associated dSij contour integral reduces, in the limit, 
to the contribution from the single pole at 6ij = 0. In such cases the prescription for setting 
some X?- — )• is therefore simply: (a) omit the factor r((5,j)P- ■ '^ from the product, (b) drop 
the integration over Sij, and (c) set 6ij = everywhere else it appears in the integrand. 

Unfortunately things are not always so simple, notably in the presence of denominator 
factors such as in the result (4.19). In that example we begin with a 9- fold integral for a 
quantity depending on 9 independent cross-ratios. Suppose we then want to take the limit 
described in the caption of Fig. 2 to recover the 4-mass double box, whose value is given 
explicitly by the formula (2.4) in terms of the 2 remaining cross-ratios which survive in this 
limit. What happened to the other 7 cross-ratios? It is easy to check that 2 of them go to 
zero in the limit. These two are simple to deal with using the rule explained in the previous 
paragraph. However, 5 of the original 9 cross-ratios go to 1 in the limit, leaving us with a 
7-fold representation for a quantity depending on only 2 cross-ratios. 

Let us conclude these remarks by noting that in examples where this kind of thing 
happens, typically it is possible to explicitly evaluate some of the 'extra' Mellin integrals 
via successive applications of Barnes lemmas. For instance in the double box case one easily 
reduce down to four Mellin integrals-the same number as the representation used in [51]. Once 
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we can no longer continue in this fashion, a possibihty would be to trade Mellin integrals for 
Euler integrals (i.e. some integrals from to 1) by using identity (B.4), since these might be 
then easier to perform explicitly. 

5 Feynman rules for conformal integrals 

5.1 Feynman rules in Mellin space 

In the last two sections we have seen that the box and double box loop integrals have Mellin 
amplitudes which are extremely simple: they are respectively 1 and 1/(1 — s), where —s is 
the square of the total Mellin momentum flowing through the internal propagator in the dual 
exchange graph {i.e. Fig 2(b)). This is in agreement with Mack's analysis of the Mellin 
amplitude, as presented in section 3. These two simple results suggest that there might exist 
some simple "Feynman rules" for writing down Mellin amplitudes. As we shortly reviewed in 
section in 3, such rules were indeed found to exist for calculations performed in AdS/CFT. 
The rules were written in a compact form in [9] and, for Mellin amplitudes of tree-level scalar 
correlators in AdS, they take the form: 

• For each internal line of conformal dimension A in the diagram write down a propagator 

(5.1) 



n\r{l + A + n-h) k'^ + A + 2n' 
In g^"^'<p'^ theory^ the vertex connecting lines with dimensions Aj is given by 






h\ 



,(m) /^_E"=l^i-2/i 



^ p(jn) I l^l^L_l , {^^,. . ., ^^^}, {l+Ai-/i,. . ., 1+A„-M, 1, • • • , 1 



with F^ one of Lauricella's multivariable hypergeometric functions. 

• The Mellin amplitude is obtained by attributing momenta to every line consistent with 
momentum conservation at every vertex and summing over all nj. 

Now consider taking the limit where the conformality condition X^I^i Aj = d is satisfied, 
keeping dimensions otherwise arbitrary. In this limit, the Lauricella function simplifies to 
one, simplifying substantially the vertex factors. However, the overall gamma function is 
divergent as well. If one goes through the computation in detail, one sees that these gamma 
functions originally arise from the integration over the radial coordinate of AdS. When the 
conformality condition is satisfied, the integral becomes divergent due to the near horizon 

®Here we really mean interactions without derivatives between arbitrary scalars e.g. (j>i(j>24'34'4, where each 
field can have a different conformal dimension. 
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region. This indicates that in this hmit, indeed the AdS integral is reducing to a boundary 
one. To deal with this divergence one can simply absorb the divergence into the definition of 
the coupling constant g^"^'. 

Since the Lauricella function has simplified, the vertex factors take a nice factorized form, 
where every line connected to it contributes a Pochhammer symbol. Therefore we may as 
well associate these factors to the propagators instead of the vertices. The modified rules 
become: 

Vertex : g'-"'^ (5.2) 

Propagator : [i, | ^ T ^^" -^ 1 (5.3) 

^^ n!r(l + A-/i) A;2 + A + 2n ^ ^ 

We see that this is almost exactly the same as in (4.16). The difference can be explained by 
different normalization conventions for propagators, and is easily ammended. With this final 
modification, we are led to the strong suspicion that the conformal flat space integrals have 
Mellin amplitudes described by Feynman-like rules, where to each internal line one associates 
a factor as in (4.16). 

Our results have direct consequences for the computation of loop integrals with dual 
conformal symmetry. They imply that for diagrams which look like tree-level position space 
correlators in cj)^ theory, one can immediately write down the corresponding Mellin amplitude. 
To reiterate the rules are: 

• To every external line in the dual diagram, attribute a Mellin momentum ki, satisfying 
kf = -1. 

• Mellin momentum is conserved at each vertex. 



To every internal line attribute a propagator factor: 

1 



p2 + l 

where p is the total Mellin momentum flowing through that line. 



(5.4) 



We should emphasize that these rules only hold for tree-level graphs in the dual position space. 
In the original momentum space variable, window-like diagrams would not be captured by 
the above set of rules. We will have more to say about this in the discussion section. Another 
important point is that these rules are really only appropriate for conformal integrals without 
nontrivial numerator factors (such as the one shown in Fig. 2.1), which so far we have not 
considered. Of course many such integrals are important for the computation of scattering 
amplitudes, and this does not mean that the Mellin space approach is useless there. Rather 
as we will see in later sections, for cases with nontrivial numerators we will end up with a 
kind of mixed representation, where the index structure of the numerators ends up separate 
from the bulk of the integral, the latter being expressed in Mellin space. 
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This is a somewhat unaesthetic state of affairs. Although we will not fully solve it in this 
note, an attempt to remedy it is given in appendix D, where we define a generalized version of 
the Mellin transform. The basic idea is to think of numerators as fields with negative (integer) 
conformal dimension. With our definition, any one-loop conformal integral with an arbitrary 
number of numerators has a generalized Mellin transform which is simply equal to one. This 
is a promising start, and we hope to explore this further in future work. For now let us turn 
to the exploration of some of the interesting consequences of the existence of Feynman rules 
for the Mellin amplitude. 

5.2 Consequences of the Feynman representation 

We have found that the Mellin amplitude that corresponds to a tree-level diagram in posi- 
tion space is given by a product of factors. For instance for diagrams in c/)'^ theory the Mellin 
amplitude becomes a product of simple propagator-like factors. This is quite nice, since prod- 
ucts of Mellin amplitudes translate into convolutions of functions in position space. A quick 
reminder of how this works in a one-dimensional example: suppose we have two functions 
f{x),g{x) with Mellin transforms M^ {s),M^{s). That is, 

/■ + 00 J /•+00 J 

Mf{s)= —x'fix), M9{s)= —x'g{x). (5.5) 

Jo X Jq X 

Then the function h{x) which corresponds to the product of the two Mellin amplitudes is 
given by 

h{x) = f^. Mf{s)M3{s)x-' = ^ ii 1^" ^ y« f{y)M3{s)x-' 

^ dy 

— f{y)9{x/y)- (5.6) 

y 

This means that we can break up the calculation of the position space functions into convo- 
lutions of simpler functions. The integral representations one obtains in this fashion can be 
thought of as solutions to certain differential equations that the original conformal integral 
satisfies. 

For clarity, consider for instance the double box integral of Fig. 2(b). We have found that 
this integral can be represented in Mellin form by 

^3,3 = i d<5., ^ n r(^«.) ^/" > (5-7) 

up to a prefactor irrelevant to the current discussion. Inspection of the diagram 2(b) shows 
that by acting with the laplacian operator on the internal line we should reduce the integral 
to something resembling a six-point star, or hexagon integral in the original momentum 
variables. In Mellin space it is clear how to see this, since polygon integrals correspond to 
setting the Mellin amplitude to one. Therefore we must come up with a differential operator 
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which cancels the propagator factor 1 — s in the above. Start by defining the homogeneous 
derivative 

^^^=P^^4^- (5-8) 

Inside the MeUin integration sign this derivative gives —Sij. Recahing (4.18) we can then 
write 

2 1 + E^*^- 43 =h = jd5^jX\T{5,,)P;^''\ (5.9) 



The conforniahy invariant functions 



Kj 



h,3iui) = (^14^25^36)43, 
hM = {PuP25P36)h (5.10) 

depend on the same cross-ratios Ui, 

Ui = Ui^i+3, i = 1,2,3, (5-11) 

Wi+3 = '"i+i,i+5, i = l,...,6 (5-12) 

with the notation of equation (2.11), 

-.. - ^^^^^^. (5.13) 

In terms of these conformally invariant functions we can write (5.9) as 

U3duj3,3 = lh- (5.14) 

Such a differential equation implies that knowing Jg should be sufficient to solve 13^3. The 
way to do this is precisely by the convolution method presented above. We shall need the 
position space expression corresponding to the Mellin amplitude 1/(1 — s): 

t<j 1^3 

with 0(x) the step function. Putting together (5.6), (5.7) and (5.15) we can write 

1 f^°° du' - 

h,3{Ui) = -- —^I(i{ui,U2,u'r^,...,Ug) (5.16) 

2 Ju3 «3 

which indeed solves (5.14). 

It is clear that this kind of strategy will generalize to the case of the L-loop ladder diagram, 
which will have L — 1 propagator factors in its Mellin amplitudes and consequently will be 
able to be expressed as an (L — l)-fold integral of the scalar 2L + 2-gon integral. It should 
also be clear that our arguments are more general, and hold not only for scalar box integrals, 
but for more complicated integrals with arbitrary legs and conformal dimensions — only the 
convolving function becomes slightly more complicated. 
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6 The chiral pentagon 

6.1 The pentagon in MeUin space 

Here we consider the chiral pentagon integral shown in Fig. 4(a), which is a basic ingredient 
in aU one- loop MHV amplitudes in SYM theory [23]. We start with the related integral 



i-^.j'i'Qi-Q-^mizhi) ''■'^ 



(the superscript indicates the presence of 1 numerator factor). Up to a certain normalization 
factor to be fixed shortly, the chiral pentagon is clearly a particular case of the above — it 
corresponds to imposing the requirements mentioned in section 2.1, 

Y-Pi = 0, i = l,...,4, (6.2) 

and further demanding that P12 = P34 = in order to match the kinematics of Fig. 4(a). 
To perform the integral (6.1), we introduce Schwinger parameters as usual, 

I'n = ^-' [fl^tt /'d^Q(-Q•y)e«■S*»^^ (6.3) 

The trick now is to trade the unexponentiated Q for a derivative with respect to the argument 
of the exponential. The calculation then proceeds in a straightforward fashion, and in fact it 
is essentially the same as for the one-loop integral without numerator which we carried out 
in section 4.1. 

With the conditions (6.2), the Mellin representation of the integral is then simply 

II = i-Y . P5) / d6,, n nS^j)Pij''' > (6-4) 

i<j 

with 

So, from the above we read off the Mellin amplitude which once again is equal to one (or, 
more precisely, (— T • P5) x 1). 

The result (6.4) is valid for general chiral pentagons. A generic five-point conformal 
integral depends on five arbitrary cross-ratios. We choose these as 

Pi,i+3 Pi+l,i+2 • _ -, r- rr.a\ 

"'^ ~ D D ' < — !,..., o y'o.o) 

^i,i+2 -fj+l,i+3 

where cyclicity is understood. If we now specialize to the kinematics of Fig. 4(a) then the only 
non-zero cross-ratios are ui^u^ and U/^. In order to compare with (2.9) it remains only to fix 
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the overall normalization of (6.1) compared to (2.7), which is uniquely fixed by normalizing 
the leading singularity to 1. This requires a little work and the nontrivial identity 

(12(345) n(678))(1247) = (1245) (1267) (3478) [1 - us - u^ + uiusm] (6.7) 

(where the momentum twistors appearing inside the brackets are labeled according to (2.7)) 
which may be derived for example with the help of appendix A of [46] . 
This leads to our final expression 

\/' = ^^14^25^35(1 -U3-Ui + UiUsUi) (f d6i, J] T{6^j) P~/'' . (6.8) 

This expression provides a simple representation of the chiral pentagon integral as a three- 
dimensional Mellin integral (from the 2^(5 — 3) = 5 original 6ij we must set S12 = S34 = 0, 
leaving only three integration variables — of course the i,j = 1,2 and i,j = 3,4 terms are 
omitted from the product). It is now trivial to verify numerically that this expression agrees 
precisely with (2.9), and in fact we derive it analytically in appendix C.3. However we will 
see in the next section that in fact the special numerator factor leads to an easier way to do 
the calculation. 

6.2 Mellin magic numerators 

Numerators of the type shown in Fig. 2.1 were called 'magic numerators' in [52] . In this section 
we will see, by careful study of the factor I — U3 — u^ + uiu^Ui appearing in the example (6.8), 
that these numerators can perform certain magic also in Mellin space. Specifically we will see 
that this factor seems carefully constructed to trivialize one of the contour integrals. After 
starting with the three-dimensional integral (6.8) and performing one of the Mellin integrals 
we obtain 






M(c3, C4) = r (1 - c-sf T (csf r (1 - c^f r (c^y 2FM, C4, 1, 1 - ui). (6.9) 

The C3 and C4 integrals run along the imaginary axis with a small positive real part. In this 
way, for [ttsl and |n4| smaller than one we can close the contour on the left and pick up the 
poles at C3 = —ni and C4 = — n2, for all ni,n2 non-negative integers. 

Now take the prefactor and place it inside the contour integral. Next perform changes 
of variables in C3 and C4 such that all terms come multiplied by uj ^^ and nj'^". Of course, 
once we do this, the contours are not necessarily the same. However, assume for a moment 
they would be. Then we are free to add up the different terms under the same integral sign, 
and we get exactly zero! This means that the original integral should be equal to the terms 
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we have lost by the shifting of the contour. This is more clearly seen with a simple example. 
Consider the integral 

.e+ioo ^^ 



il -.>[•_ 



2TTi 



r{sfr{i-sfx- 



(6.10) 

' t—ioo -^^^ 

with < e < 1. For < |x| < 1, we close the contour on the left and pick up the poles at 
s = —n, with n a non-negative integer. The result is found to be — log(x). Now let us do the 
calculation in a different way. By placing the prefactor inside the integral and doing a change 
of variables we find 



(1-x) / 



e—ioo 
e+ioo 



I^^<"^^(1 



s? X-' 



/ f.Tisfrii-sfx-^ + 

Je-ioo -^^^ Je-1 



e— 1+ioo 



(/■£— 1+ioo 
/ 
J e — l — ico 

-log(x). 



ZTTl 



e+ioo 



d^ 

2Td 



r(s)^r(i-s)^ X' 



(6.11) 



In the last step we used the fact that the contour in the last integral only encloses the pole 
at s = 0. In other words, in this case the prefactor is so special that the result of the integral 
becomes the residue of a single pole, as opposed to an infinite sum. 

This is precisely what is going on in our more complicated example of the pentagon 
integral (6.8). There we have four terms coming from the prefactor 1 — its — M4 + U1U3U4. 
Multiplication by the last three of these four terms involves shifting the contour, and just as 
in our simple example we get the zero mentioned previously plus the terms coming from the 
contour shift. The result is 



y27^ = |l^esc3=o / ^.u^'^'u^'^Mica + 1,04) 



-Ul ReSc4=0 / ^^3 ^«U4 ^*M(C3 + 1, C4 + 1) + (C3 o C4) !> + (6.12) 

+Mi ReSc3=o,c4=o u^''^u^'''^M{c3 + 1, C4 + 1). 



(6.13) 



Evaluating the residues leads to 



log('Ui)log(u3)+log(u3)log(n4)+log(u4)log(ni) + Ul 



d^ 



2Fi(a,/3, 1,1-ni) 



/^(%^+V)r(i-c)2r(c) 



2 A 

da 



da dp 
2Fi{c,a, 1, l-ui)-ui2Fi{l+c,a, 1, 1-ui] 



a=l3=l 

(6.14) 
0=1 



To proceed we need to evaluate the derivatives of the hypergeometric function with respect to 
its parameters. Such derivatives can be easily evaluated by first using the series representation 
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for the hypergeometric function, taking the derivative, and performing the sum. In this way 
we find 



da 



2Fi{c,a, 1, 1 - Ml) - Ml 2-^1(1 + c, a, 1, 1 - mi) 



1-M^'= 



a=l 



2Fi{a, (3,1,1 -Ml) 



= ^(1-Mi)"i/^ 
a=l3=l n=0 



log(Mi)2 + Li2(l-Mi) 



Ml 



(6.15) 
(6.16) 



where Hn = X]fc=i ^ ^ is the n-th harmonic number. We get 
dc 



2iTi 



(^- - 1)(^3- + u^-) r (1 - c) r (cY r (-c) 



+ log(Ml) log(M3) + log(M3) log(M4) + log(M4) log(Ml) + log(Ml)^ + Li2(l - Ml). (6.17) 

The single Melhn integral can be easily performed and we recover the well-known result for 
Is shown in (2.9). 

To conclude this section, we recall that in [50] it was found that the pentagon integral 
satisfies the particularly simple differential equation 



M3M4a„35„,/« 



1. 



(6.18) 



Using the representation (6.17) it is trivial to see this. Clearly this operator annihilates the 
single Mellin integral, since in it M3 and M4 appeared summed, not multiplied. The action of 
the operator on the remaining terms immediately gives 1. We note that to obtain this result 
it is actually not necessary to compute the derivatives of hypergeometric functions. One can 
start directly from expression (6.9), apply the differential operator, and shift contours. One 
gets an expression analogous to (6.13), and quite trivially one immediately obtains one as the 
answer. 

7 The chiral hexagon and differential operators 

Now we consider the hexagon in Fig. 4(b). In the ambient space formalism this integral 
translates into (again, up to an overall normalization factor which will be fixed below) 



where the special vectors Y, Y' satisfy 



^^K-Q-Y){-Q-Y') 



d^Q 



nUi-p^-Q) 



Y-Pi = 0, i = l,...,4, 
Y' -Pi = 0, i = 3, . . . , 6, 



(7.1) 



(7.2) 
(7.3) 



and further we demand -Pi,j+i = for all i. There are two possible solutions to each of the 
above equations, for a total of four different integrals. However two of these are related to the 
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other two, leaving only two independent choices. Our focus here will be on the chiral integral 
shown in in Fig 4(b) and (2.8) in which both numerators are of the same type. 

As always we begin with a more general n-gon integral, this time with two numerators, 



r 



yAyB 



vr" 



^^QQaQbW 



r(A, 



\[-Pi-Q)-^^' 



(7.4) 



and again we assume the conformality condition ^" Aj = d+2, but keep all vectors arbitrary. 
To do this integral we follow the same strategy as for the pentagon, obtaining: 



with 



i2 = y^y^/d5.,r^^nr(^^^)^7 



TAB ^ ^AB + Y,5.^^ + Y.P^^PFS.. 



(7.5) 



i¥^3 



(7.6) 



«j 



Here we have defined Si as the operator which modifies the constraints on the 5ij in such a 
way that Aj — )• Aj + 2, keeping all other dimensions fixed. This somewhat technical point 
will not be relevant for the chiral hexagon. Notice that Tab is symmetric and traceless as it 
should be, due to the fact that 



^AB 



tj^^'^Tab 



d + 2-^5ij 



d + 2 - ^ Aj 



0. 



(7.7) 



j=i 



There is a very nice way of rewriting equation (7.5). To see this notice that inside the Mellin 



integral, we have 6i 



P.. d 



-din. Therefore we can write 



ll=YAY^\r^ 



AB 



pApB 



i,j 



Pi. 



^^,+Y,P^^PFS^ 



<^^^J\{n5^,)P^ 






i<j 



(7.8) 



That is, the two-numerator polygon conformal integral can be written as a certain differential 
operator acting on the expression between square brackets. The object being acted upon 
looks exactly like the polygon integral without numerators, namely (4.10). The only catch is 
that for the polygon integral we have ^" Aj = d, whereas here, on account of the original two 
numerators, we must have ^^ Aj = d + 2 for conformality. Therefore the object appearing 
above could never arise from the computation of a d-dimensional conformal integral. Instead 
it must come from a. d + 2-dimensional integral dependent on the same collection of cross- 
ratios (which are largely ignorant of the dimensionality of the space they live in)^*^. This is 



^"They may satisfy polynomial relations known as Gram determinant constraints, but these have no effect 
on any of our analysis. 
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the Mellin space manifestation of the weh-known connection between scalar integrals in d + 2 
dimensions and tensor integrals in d dimensions (see for example [60]). 

Let us now go back to the particular case of the chiral hexagon in d = 4. The result 
can be written as a particular operator acting on the hexagon integral in d = 6. The generic 
integral depends on nine independent conformal cross-ratios, as introduced in equations (5.11) 
and (5.12). In the kinematic regime relevant to the chiral hexagon we have Pi.i+i = for 
all i, and so the cross-ratios M4 through ug actually vanish. The Mellin representation of the 
integral (2.8) including all overall factors, and with the numerator factor appearing in (7.5) 
and (7.6) written out explicitly in terms of momentum twistors is given by 

3 4 




P15P24P36 i> d<5i,- M{5,,)J{T{5i_ 

i<j 






(7.9) 



1 6 

where (see appendix C.4 for details) 

(1345) (2346) ^ (1356) (2346) ^ (1246) (1345) ^ (1246) (1356) ^ 
^^^'^^ = ^- (1346) (2345) '^^^- (1346) (2356) '^^^- (1245) (1346) '^^^- (1256) (1346) '^^^ ^^'^^^ 

and of the 2^(6 — 3) = 9, the six of the form 5j_i+i are set to zero (and the corresponding 
terms omitted from the product), leaving only three independent integration variables. We 
have checked numerically that the expression (7.9) agrees precisely with (2.10). 

Alternatively, we can also derive the differential operator form of the chiral hexagon. 
Since the integral (7.1) is conformally invariant, we can trade derivatives with respect to Ptj 
by derivatives with respect to cross-ratios, and in the end we obtain 



U2{l-Ul-U^)+U2{l-Ul)uidu^ 




1 6 

+U2{l-u-i)u^dus-{l-U2){l-ui-u-i)u2duATQ{ui,U2,uz) (7.11) 

where the zero-mass d = 6 hexagon 

Mui,U2,u^) = {Pi^P2^Pmh)u,=...=u,=^ (7.12) 

has been given explicitly^^ in [38] (but see also appendix C.2). Of course, this particular 
differential equation relating the chiral two-numerator hexagon in Ad to the massless scalar 
hexagon in 6d is very well-known [50] (and see also [37] for further applications). 



^'^ Regrettably, the notations of [38] and [23] are inconsistent. We stick with the latter, so Del Duca et al's 
Mi,«2,W3 are our U2,Ui,u\ respectively. 
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Figure 5. The simplest example of a dual conformal 'window' diagram (black) whose dual diagram 
(blue) has a loop in x-space. Here we have shown the fully massive version of the integral; the fully 
massless version contributes to the four-loop four-particle MHV amplitude in SYM theory. 

8 Conclusion and discussion 

Motivated by the success of Mellin representations [7] for studying correlation functions in 
general CFTs and in particular in AdS/CFT, we have here initiated a preliminary investi- 
gation of the suitability of using Mellin representations for dual conformal integrals of the 
type appearing in SYM theory scattering amplitudes'^. We have explored this possibility by 
working out explicit, and very simple, Mellin amplitudes for several particular integrals. In 
those cases for which results are available in the literature, including the four-mass box, the 
four-mass double box, and the chiral pentagon and hexagons, we have checked agreement 
(either analytically or numerically) between these results and our Mellin representations. Be- 
yond this, we have seen how simple it is to write Mellin representations for large classes of 
integrals which seem far beyond the ability to evaluate in u-space with currently available 
methods. Examples of this include the fully massive L-loop ladder diagram, or even the 
fully massive one-loop n-gon integral in n dimensions'^, whose Mellin representation, when 
normalized according to the standard convention (3.5), is just l! 

This was made possible by introducing a set of Feynman rules for conformal integrals, 
and showing how they can be interpreted as arising from the Feynman rules found in [6, 9] 
for AdS correlation functions. Unfortunately, this analogy with momentum space amplitudes 
does not hold at the quantum level, which is to say whenever there are loops in x-space 
(window- like diagrams in the original momentum space, see Fig. 8). This is because as we 
saw in section 3, the Mellin amplitude can have only simple poles, whereas momentum space 
loops have branch cuts. Of course this doesn't mean that there aren't some simple rules for 
writing down Mellin amplitudes involving x-space loops — indeed it is extremely likely those 
must exist — it is just that they have yet to be worked out in detail. For some work along this 
direction in the context of AdS computations see [8, 16]. 



^■^As opposed to more general Mellin-Barnes representations, which are a general tool for a much wider class 
of Feynman integrals. 

^^For even n, an explicit formula for the symbol of this integral is given in [40], however the fully massive 
version of the integral has not yet been evaluated for any n > 4. 



27 



It would be particularly interesting to work out if there exist similar diagrammatic rules 
for integrals involving general numerators of the type reviewed in section 2.1. In the examples 
we have studied we have seen that for integrals with only a single numerator, the numerator of 
the Mellin representation is just a constant prefactor (independent of the integration variables 
5ij), while for integrals with more than one numerator factor a non-trivial numerator can 
appear such as in (7.9). Because these numerator factors depend explicitly on cross-ratios, 
these Mellin amplitudes are not literally the same thing one would obtain if one computed the 
inverse Mellin transform of the integral's u-space answer. Rather they are some kind of hybrid 
representation (akin to writing a function of x as a Fourier transform of something which 
depends on k but also explicitly on x). This is not necessarily a bad thing; if for example there 
are simple and physically well-motivated general rules for writing such hybrid representations, 
then we are all for it. Alternatively, the generalized Mellin transform introduced in appendix D 
might also be a natural object to consider. 

A particularly important role is played by the one-loop scalar n-dimensional n-gon in- 
tegral, whose Mellin amplitude is exactly 1 (in the normalization which is now standard in 
the AdS /CFT literature). We have seen in the concrete example of the chiral hexagon how 
the numerator factors transform into certain differential operators acting on these simpler 
integrals. Generically, the statement is that a d-dimensional conformal integral containing an 
even number 2n of numerators can always be rewritten as an operator acting on the same 
conformal integral without numerators in a higher dimension. This operator contains up to 
n derivatives with respect to the Pij, and in turn these will transform into derivatives with 
respect to cross-ratios. In this way Mellin space allows for an easy deduction of differential 
relations amongst multi-loop integrals. 

Going in the other direction, we have shown that certain higher-loop integrals with trivial 
numerator factors may be written as integral operators acting on the n-gon. Our chief example 
is a formula for the 6-mass two- loop double box integral as a very simple integral of the 6-mass 
scalar hexagon in six dimensions, 



."•■■■'"-^/rt^^Yiy^ *"'■■■' •'■'' 



where u = xI^xIq/^xI^x^q), found in section 5.2. Here the diagrams on the left- and right- 
hand sides are each functions of 9 cross-ratios, eight of which (represented by . . .) are identified 
on both sides. It is only the appearance of one preferred u which breaks the cyclic symmetry 
of the hexagon integral. Again we emphasize that explicit formulas in u-space are currently 
not known for either diagram, but the double box integral on the left does appear in two-loop 
SYM theory scattering amplitudes. This formula clearly generalizes to one which expresses 
the fully massive L-loop ladder diagram as an L — 1-fold integral of the fully massive scalar 
2L + 2-gon. 
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Another natural arena in which to explore these methods is in the context of the chiral 
double pentagon integral, which was evaluated in [44] as an ingredient in the analytic for- 
mula for the two-loop 6-point NMHV amplitude presented in that paper. It satisfies known 
differential equations relating it to both the chiral hexagon and the six-dimensional scalar 
hexagon [37, 50], which manifest themselves very simply in Mellin space. Based on our anal- 
ysis we expect that it should also be possible to express the double pentagon as a first-order 
differential operator acting on a single integral of the 8-dimensional octagon. This relation is 
especially interesting in light of the special role played by a set of chiral octagon integrals [23] 
which provide a basis for one-loop integrands and have been evaluated analytically. These 
chiral octagons should themselves be expressible as certain differential operators acting on the 
8-dimensional octagon. This suggests a possible link between the double pentagon integral 
and these chiral octagons, which would surely be worth further investigagion. 

To summarize, it is clearly important to better understand the structure of the scalar 
n-gon integral in n dimensions, whose Mellin amplitude is just 1. Since all numerator factors 
in a general integral can be written as differential operators acting on this object, and all 
denominator factors can be written as integral operators acting on this object, we believe 
that all (at least, all non-window) dual conformal integrals relevant to SYM theory can be 
written as integro-differential operators acting on the scalar n-dimensional n-gon. 

On a different note, it would be very interesting to understand, for those integrals which 
can be expressed in terms of generalized polylogarithm functions, if there is a simple way to 
read off the symbol of an integral directly from its Mellin representation (or even to be able 
to look at a Mellin representation and quickly determine what degree of transcendent ality, 
if any, the resulting integral has). Presumably, if there is a way to determine this it will be 
by careful examination of the arguments of the gamma functions which appear generically 
in any Mellin transform. Indeed, (poly)logarithms can only arise from sums of double or 
higher order poles in the Mellin amplitude. A first step would surely be to understand the 
n-dimensional scalar n-gon integrals, which have Mellin transform 1 but a very nontrivial 
symbol [40]. 

Keeping our wildest speculation for the very end, it is perhaps our greatest hope that 
it might be possible to find a recursion relation which works directly at the level of Mellin 
representations for multi-loop amplitudes in SYM theory, akin to the BCFW-type recursion 
relation which is known to hold at the level of the integrand [28-30]. If this hope was 
realized, and one could generate Mellin representations for arbitrary amplitudes 'at the touch 
of a button', then our suggestion that Mellin space might serve as a useful stepping stone 
between integrands and integrals would be fully realized. 

A The Symanzik star formula 

For completeness, in this section we review the Symanzik star integration formula in Euclidean 
space as discussed in [7]. For a proof and more details we refer the reader to the original 
reference [58] . Consider a set of n points in Euclidean space Xi and their differences Xi — Xj . 
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In the embedding formalism we have Pij = —Pi-Pj = (xj — Xj)"^. Symanzik's formula is then 

P(nfi^)e-(s^.-.."*'*^^-^) = ifd5., n m,)/^7^-. (A.i) 

This identity can be deduced by using the Cahen-Mellin integral e~^ = J i^ r(s) x~^ for 
each factor exp{titjQij) and performing the Schwinger parameter integrals. The integration 
measure on the right-hand side is more precisely given by 

/p+ioo J^ 

dSij = 2 H d5ij n SiAi - Y, Sij). (A.2) 

We have included the factor of two since in the majority of cases one solves for a subset of the 
6ij as independent variables, and in that case it is easy to see that a factor of 1/2 arises from 
the Jacobian coming from the delta functions. Of course, more general choices are possible. 
To see this, first notice that the Dirac delta functions imply constraints on the 5ij: 

Y^S^^ = Aj = -5ii (A.3) 

for all i. Now pick a particular solution of the set of equations (A.3), 6^,. Then we write 

|n(n— 3) 

Sij = <5°- + Yl ^v,kSk (A.4) 

with 

Cii,k = 0, Y ^ij,k = 0. (A. 5) 

Choosing as independent variables the (yi2n{n — 3)) coefficients Cij^k with 2 < i < j < n (with 
the exception of 023,^)1 with the further restriction |detcjj_fc| = 1, we can write 

, . in(n— 3) 

/+JOO 2 >. ' He, 
.. n S- (*■"' 

"*°^ fc=i 

The integration paths are chosen parallel to the imaginary axis, with real parts such that the 
real parts of the arguments of the gamma functions are positive. 
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B Useful hypergeometric function identities 

In the appendices we show how one can (sometimes) easily go from the Melhn representation 
to position space. We will be using the following basic identities 

r(c) /" d« .„ ^,,_s, 



^^'("•"■^'^'^r^SjfiFT^X ^j^f(i-*r»(i-*.r, (B,i) 

^, , , f ds T(c)T(s)T(c — a — b + s)T(a — s)r(b — s) , , „ ,^ x 

.F.(a, i.c..)^f- 'rl^rimc - 4r(c - /) " ' ^' <''■'' 

r. / , X r(c) / ds r(s)r(a - s)r(b - s) , , , ,^ , 

From the first two equations above we also get the useful relation 
/ ^ r(ai + s)T{a2 + s)r(6i - s)r(62 - s) z"^ 

= r(ai + 6i)r(a2 + 62) / ^7^^ t'^^'^^ (1 - t)"^+'^ [1 - *(! - z)]-''^-"^ . (B.4) 

Jo f-{^ ~ <^) 

which allows to go from a Mellin-type integral to an Euler-type one. 

C Selected details 

C.l The general box integral 

The box integral is written in Mellin form as 

4 



/4=y"d<5., nm,)/^/-. (C.l) 



Defining the cross-ratios 
and the quantities 
the integral becomes 



Kj 



P12 PsA Pu P23 ,r. oA 

A3 ^4 -ri3 ^24 

AJ.^. A..^. ,C.3, 



with 

^^r(A+ + s)r(A+ + s)T{-A^^ - s - t)r(A3-4 -s-t) 

X r(t)r(A^2 - ^34 + 1) u~' v~\ 
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Using the integral representations for the hypergeometric function it is straightforward to 
obtain 

f(^^) = r(A+ -Ar^) r(A+ -A3-J r(A+ +Ar2) r(A3+ +A3-J r' _^^ a_t)^t.-^i2 
^^ ' ' r(A+ +A3+J Jo t{i-t) ^ ' 

X tA^.+A3- (,_,(,_^))-A^,+A- ,i.,(^A+-Ar2,A3++A3-4,A++A+,l-i^^). 

The above integral can be done in terms of Appell -F4 functions, but we will not do so here. 
Instead, we specialize to the case of equal conformal dimensions, A,j = n. Then we obtain 

.. , r(n)4 /■! dt / t(i-t) Y r. ( {l-t)tv 

fA{u,v)^-^-^ — \.^ ' , 2i^i n,n,2n,l ^ ' 



r{2n) Jo t{l-t) \l-t{l-ui)J \ ' ' ' l-t{l-u) 

r(n)4 /■! dt -n r, f 01 l-til-u)\ ,^^, 

V '^ 2Fi{n,n,2n,l ,/ . ■ (C.5) 



r(2n)yo t{l-t) ^ ^ \^ ' ' "' t{l-t)v 

The case relevant for d = 4 is n = 1, whereupon the integral reduces to 

It is now a straightforward matter to integrate this expression to obtain the final result (2.4). 

C.2 The six-dimensional one- mass hexagon 

In this section we derive a simple Mellin-Barnes representation for the hexagon integral in 6 
dimensions. The Mellin amplitude is one, and so the integral which we denote Iq is given by 

h = fd5,jl[T{6,,)p-/^^ (C.7) 

i<j 

By conformality, up to some prefactor, the six-point integral can depend only on 6 x 3/2 = 9 
cross-ratios Ui. In the notation of (2.11) we take these to be given by 

Ui = Ui^i+3, « = 1,2,3 (C.8) 

Ui+3 = Ui+i,i+5, i = l,...,6. (C.9) 

Accordingly the constraints on 5ij can be solved for nine independent variables which we 
denote by Cj. There are many possible choices for which variables to choose, and this cor- 
responds to shifting the function of cross-ratios by some monomial in the cross-ratios. One 
particular choice leads to 

h = {PUP25P36) h= flip: u;"^ ( n r(Q) ) r(ci - cs - ce)r{c2 - ce - C7) 

Xr(l - Cl - C3 + C5 + Cs) r(c3 - C7 - Cs) T{ci - Cg - Cg) r(l - Cl - C2 + Cq + Cg) 
Xr(l -C2-C3 + C4 + C7) r(c3 - C4 - C5) r(c2 - C4 - Cg), (CIO) 
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To proceed we will focus on a specific kinematic regime, corresponding to the one-mass case. 
In this limit all cross-ratios except ui, U2, U3 and U4 vanish. The corresponding Mellin integrals 
become trivial, as one simply picks up one residue for each of the cross-ratios. The integral 
becomes 

^6 = - / (n ^ n-^" J r(ci)2 r(c2)r(c3)r(c4) r(i - ci - 02) 

X r(l - Cl - C3) r(c2 - C4) r(c3 - C4)r(l - C2 - C3 -F C4). 

By making good use of the identities of section B we can simplify this to 

dC2dC3 _e2 _,3r(l-C2)3r(c2)2r(l-C3)3r(c3)2 



(2vri)2 "2 "3 r^2 - C2 - C3) 

X2F1 (1-C2,l-C3,2-C2-C3,l-Ui) 2F1 (C2,C3,1,1 - Ui) . 

This provides an efficient numerical representation for the integral. Of course nothing stops 
us from trading the Mellin-Barnes integrals for Euler integrals. This can be done for instance 
by using the Euler representation of the two hypergeometric functions and performing the 
Mellin-Barnes integrals. This leads to the alternative representation 



^ ^°g [ {l-s){l-t{l-U4))u2 



l-s(l--ui) 
dsdt- ,.. .)il-til-u,-, __ ^^_^^^ 

° ' SUl + (1 - S)[l - (1 - t(l - U4))U2]]I{1 - t)sU3 + t 



We have checked that this agrees numerically with the exact analytical result of reference [38] , 
which was denoted there by Xg^^ In doing the comparison note that our cross-ratios correspond 
to the ones in that reference after we change variables such that 

h{ui,U2,Uz,Ui) =X6,m(^2,'"3,Ui,U4). (C.12) 

C.3 The chiral pentagon 

There are five Mellin integrals in the expression for I^. After solving the constraints on the 
6ij we find we can write 



-P25 -P14 -P35 U ■ 1 

2 = 1 



/ n (^<'"'r(l - ^^)) r ^^1 + C2 - C4 - 1) 

X r (C3 -h C4 - Cl) r (C2 + C3 - C5) r (Ci -h Cs - C3 - l) r (C4 -h C5 - C2) 



with the cross-ratios 



Ui = Ui^i+2 = -p, 5 > ^ = l,...,5. (C.13) 
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We are interested in evaluating the Mellin integrals for the particular kinematics correspond- 
ing to i-*i2 = -P34 = (as in Fig. 4(a)), which implies U2 = u^ = 0. In this limit we simply 
pick up the residues at C2 = C5 = 1, obtaining 

/ dcidc3dc4 I I c4-ir/i \-r(-, \-r (-, ^ 

/ (27ri)3 ^ 3 ^ r(l-ci)r(l-C3)r(l-C4) 

r (ci - C3) r (C3) r (ci - c^) r (04) r (c3 + C4 - ci) . 

Now we perform the ci integral using the identity (B.2), to get 

/ tS^'"''''''^'''^ ^^ ~ '''^' ^ ^'''^' ^ ^^ ~ '''^' ^ ^'''^' '^'^^ - C3, 1 - C4, 1, 1 - Ul). 

By further using the representation (B.l) for the hypergeometric function, the C3 and C4 
integrals can also be performed, finally giving 

^^^ log(u3[l-t(l-ni)]) 

{l-U-i[l-t{l-Ui)]){Ui{l-t) + t) 



1 



1 — U3 — li4 + 141^3^4 



1 lil 

log(ui) log ( ) - log(n3) log(ti4) - Li2(l - U3) 

1 — UiUi 



+ Li2(l-n,.3)-Li2 (^^^] +Li2 (--^^^ 



(C.14) 



One can check numerically that the expression between brackets is exactly (2.9), and that 
the prefactor in front of the expression in square brackets is exactly the appropriate relative 
normalization factor appearing as the overall factor in (6.8). 

C.4 The Mellin numerator of the chiral hexagon 

Here we provide some details on the derivation of the Mellin numerator (7.10) for the chiral 
hexagon integral. We begin by solving the constraints (7.2), which in x-coordinates read 

(y-Xi)2 = 0, i = 1,2,3,4, (C.15) 

{y'-Xif = ^, i = 3,4,5,6. (C.16) 

Details on solving this kind of problem in general may be found in [23]. There are two 
solutions for each of y and y': 

y = (13) or (612) n (234), (C.17) 

y' = (46) or (345) n (561), (C.18) 

where we express the 4-dimensional vectors as antisymmetric products of momentum twistors 
as in (2.16), labeling the legs and faces of the hexagon as shown in (2.8) and Fig. 4(b) 
respectively. Amongst the four choices of {y, y'} two pairs are related by parity. The chiral 
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hexagon we are interested in corresponds to choosing the same type of solution for each 
numerator, so we wih proceed with 

y = (612) n (234), y' = (345) n (561). (C.19) 

Now we turn our attention to the relevant portion of the prefactor in (7.8), which we can 
write as 



/\2 



(y-y) 



^ {y-y'y{xi -XjY *^ 



(C.20) 



Only four of the cross-ratios appearing in the sum are nonzero, and those are easily computed 
using (C.19), which gives exactly (7.10) together with the overall factor 

(y - y'f = (1256) (1346) (2345) = {-Pi^){+Pm){-P2i) (C.21) 

appearing in (7.9). It may be of interest to note that two of these cross-ratios can be written 
easily in terms of the standard Ui defined in (5.11) and (2.11): 

(1356) (2346) _ _ (1246) (1345) _ 

(1346) (2356) " '' (1245) (1346) " '' ^^ > 

The other two cross-ratios (the first and fourth in (7.10)) are the two roots of the quadratic 
equation 

U2X^ - {I-U1 + U2- u^)x + {I- ui-uz + U1U3) = 0. (C.23) 

D A generalized Mellin transform 

Suppose one wants to compute an integral with a general numerator structure of the form 

tN,M ^ Y[i=l^i^i) f ,d^ Y[a=liQ • ^g)'"" /-p, , X 

The trick is to introduce Schwinger parameters for the denominators and complex integrals 
for numerators, 

I-'- = i f^'Q £^Ilftfl[ ^^^a'-^ exp[Q.(j:t.P. + 5:..y.)]. (D.2) 



i=l 



The integral in Q is then simply performed. Conformality makes annoying factors of ^ tj -|- 
^ Zj drop out. The net result is that we get 

,N,M = 2£^ii ftt / n (g^ .7-- exp[(^ UP. + Y. ^.y.f\ ■ (D-3) 



i=l 



Now we do Symanzik's trick, slightly generalized. We slice up the Pij exponentials and 
introduce Mellin integrations 5ij for each of them, as in A. For the cross terms tiZaPi ■ Ya and 
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for ZaZfjYaYiy we substitute the exponentials by their series representation, introducing sums 
over integers qia and n^fe. Notice that the parameter qij has indices running over different 
ranges. The integration over the t's gives Dirac deltas imposing constraints on 6ij Mellin 
parameters; the complex Zi integrations give Kronecker deltas imposing constraints on the 
Uab parameters and we get 

TV M 

i^M = Y, f ^^^J n r(^^^) (^^ • p^y'^' n ^ (^« • ^f*)""' n ^ (^* • ^-)'"" • (o-^) 

„ „ J ■ , ■ ,, ^ab- . Qia- 

nabyQia i<3 a<b i,a 

The sums and the integral satisfy the constraints 

j¥=i <* 

^"-afe + ^gia ='m-a. (D.6) 

bf^a i 

This example suggests it is natural to introduce a generalized Mellin transform for any con- 
formal integral, by simply adding a general function M{5ij,qia,nab) in the integrand above. 
The generalized Mellin amplitude of any of the conformal integrals /^'^ is one. With similar 
calculations it is not too hard to show that the Mellin amplitude of an exchange integral is a 
simple pole, but the location of the pole depends on the quantities Uab and qia- In other words, 
everything works out as if we could attribute negative conformal dimensions to numerator 
factors. 
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